Quantum super-cavity with atomic mirrors 
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We study single-photon transport in an array of coupled microcavities where two two-level atomic 
systems are embedded in two separate cavities of the array. We find that a single-photon can be 
totally reflected by a single two-level system. However, two separate two-level systems can also 
create, between them, single-photon quasi-bound states. Therefore, a single two-level system in the 
cavity array can act as a mirror while a different type of cavity can be formed by using two two-level 
systems, acting as tunable "mirrors", inside two separate cavities in the array. In analogy with 
superlattices in solid state, we call this new "cavity inside a coupled-cavity array" a super-cavity. 
This supercavity is the quantum analog of Fabry-Perot interferometers. Moreover, we show that 
the physical properties of this quantum super-cavity can be adjusted by changing the frequencies of 
these two-level systems. 

PACS numbers: 32.80.Qk, 73.22.Dj, 42.50.Pq, 85.85.-(-j 
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I. INTRODUCTION 

In quantum networks, photons provide faithful quan- 
tum information transfer, because they travel at the 
speed of light over long distances, and with little decoher- 
ence compared to other information carriers (e.g., elec- 
trons). To interconnect networks, it is crucial to have 
a quantum memory at the switching nodes. Many ap- 
proaches have been proposed to realize quantum memo- 
ries, where quantum information can be stored and re- 
trieved, for instance, using electromagnetically-induced 
transparencyfe^., in Refs. [U, H, Q) or photon echoes 
(e.g., Refs. d, H, Photons can be confined to a 

very small volume (e.g., Ref. 0]) using micro-cavities or 
micro-resonators with low dissipation and thus the micro- 
cavities can serve as quantum memories. Moreover, ex- 
periments also demonstrated that the quality factor of 
a photonic crystal nanocavity [1, [§] or microwave cav- 
ity [U can be controlled by dynamically changing 
the environment of the cavity. 

To faithfully transfer quantum information, individ- 
ual photon control would be desirable. Single-photon 
turnstiles have been studied in, e.g., Refs. [ll, fisl . [l^ . 
There, a semiconductor quantum dot 12] or a sin- 
gle atom [ij, [13] can behave as a photon turnstile. 
Recently, a nonlinear two-photon switch device using 
nanoscale surface plasmons has been [l^ theoretically 
studied, where a single-photon 'gate' is used to control 
the propagation of subsequent single photons. Consid- 
ering the one-dimensional scattering process of single- 
photons by a two-level system, the total reflection can 
be controlled by tuning the inner structure of the scat- 
terer [H, [13, [ll|. A solid state device, functioning 
as a single-photon quantum switch [l^, in a one- 
dimensional coupled-cavity waveguide has been studied, 
also using a discrete-coordinate approach [l7| . 



For single-photon transport in a one-dimensional 
waveguide, the photons can be totally reflected [H, [l3, 
[Tsl by a controllable two-level system which can act as a 
perfect mirror. It is known that the Fabry-Perot cavity, 
which consists of two highly reflecting planar mirrors, is 
the simplest cavity. It is then natural to ask the ques- 
tion: "is it possible to construct a quantum resonator, in 
a one- dimensional waveguide, with two controllable quan- 
tum scatterers?^^ Here, we focus on this question and 
study quantum analogs of the Fabry-Perot cavity. 

This paper studies the coherent transport of photons, 
which propagate in a one-dimensional coupled-resonator 
waveguide (CRW) and are scattered by two controllable 
two-level systems located separately in the CRW. Be- 
sides presenting a unified theory, including both the long- 
wavelength and short-wavelength regimes, the discrete 
coordinate approach employed in this work shows that: 

1. Photon quasi-bound states, with a tunable leakage, 
appear in the region sandwiched between the two 
two-level systems, when the interaction between 
the two-level systems and the cavity field is strong 
compared with the hopping constant. (Hereafter, 
for brevity, we will often use the word "atom" in- 
stead of "two-level system". In the terminology 
of quantum information, a two-level system is a 
qubit, therefore, hereafter, a two-level system is 
sometimes called a qubit, but this "atom" refers 
to an artificial atom made, e.g., from a supercon- 
ducting circuit.) 

2. A perfect quantum super-cavity, confining photons 
inside the two "atoms", can be formed when the 
transition energies of these two-level systems are 
equal to the photon energy, with wave number k — 
n7r/(2d), where 2d is the distance between the two 
atomic scatterers and n is an integer. 



2 




The CRW can be described by the Hamiltonian 



FIG. 1: (Color online). Schematic configuration for a quan- 
tum super-cavity realized by two atoms embedded in two sep- 
arated cavities of a coupled resonator waveguide, as shown in 
(a). Each atom can be represented as in (b). These two atoms 
can behave as two partly-reflecting mirrors, forming a cavity- 
within-a-cavity, or super-cavity. This can act as a quantum 
analog of Fabry-Perot interferometer. 



3. Photons can be stored and re-emitted by adjust- 
ing the transition frequencies of these two two-level 
systems. 

This paper is organized as follows: in Sec. |TT1 we 
present our model, a coupled-cavity array with two atoms 
separately embedded in two different cavities. In Sec. lIIIi 
we study the transport properties of a single-photon, and 
derive the conditions for the coherent control of a single- 
photon scattering by two atoms. In Sec. IIVI the quan- 
tum super-cavity, with two atomic mirrors, is studied. 
We prove that the leakage of this super-cavity is tunable 
by changing the transition energy of these two atoms. 
The wave numbers inside this super-cavity are also an- 
alytically obtained by a perturbation approach. More- 
over, we study how a super-cavity can be formed in the 
long-wavelength (low-energy) regime, in Sec. |V] and in 
the short-wavelength (higher-energy) regime, in Sec. IVII 
These two regimes correspond to the quadratic and lin- 
ear photon dispersion relation, respectively. Conclusions 
are summarized at the end. 



II. MODEL 

As shown in Fig. [U we consider a one-dimensional 
coupled-resonator waveguide (CRW) with two two-level 
systems, embedded separately in two distant cavities. 
The CRW can be either an array of coupled super- 
conducting transmission line resonators or an array of 
coupled defect resonators in photonic crystals (see, e.g. 
Ref. [2l|). The two- level systems can be either natural 
atoms or artificial atoms (e.g., superconducting qubits or 
semiconducting quantum dots). In contrast to the simi- 
lar configurations in Refs. [H [H, [H, IH, , here only 
two atoms are located inside the CRW. 

Once a photon is inside one cavity of the CRW, it prop- 
agates along the CRW and is also scattered by the atoms. 



j = -oo 



j=-oo 



h.c. 



(1) 



with the annihilation operator aj and creation operator 
a'j of the jth cavity mode. The first term of Eq. ([1]) 
denotes the free Hamiltonian of all the resonators. The 
second term of Eq. ^ represents the couplings between 
any two nearest-neighbor cavities. For example, aja]_|_]^ 
means that the photon is annihilated in the jth cavity 
and is created in the {j + l)th cavity. Here, for simplicity, 
we assume that all resonators have the same frequency 
uj and the hopping energies ^ between any two nearest- 
neighbor cavities are the same. The hopping energy ^ is 
determined by the inter-cavity coupling. He is a typical 
tight-binding boson model and can be rewritten as 



bibk 



(2) 



by introducing the Fourier transform 



(3) 



The dispersion relation between and k is given by 



Ek = oj — cos fc. 



(4) 



which forms an energy band. Here, the lattice constant 
I is set to unity. 

Let us assume that each atom has a ground state \g) 
and an excited state |e). Let the distance between the 
two atoms embedded in the CRW be 2d. For convenience, 
we take the 0th cavity as the coordinate-axis origin. As 
shown in Fig. [1] we also assume that the first atom, with 
transition energy fii, is located at the (— c?)th cavity, on 
the left side of the origin, and the second atom, with tran- 
sition energy 172, is embedded in the dth cavity, on the 
right side of the origin. Under the rotating-wave approx- 
imation, the interaction between the dth and (— d)th cav- 
ity fields and the two atoms is described by the Jaynes- 
Cummings Hamiltonian 



Hr 



E 

1 = 1,2 



Qi |e), {e\ + Ji (\e)i {g\a^_j^y^ + 



h.c 



(5) 



where J; is the coupling strength between the ^th atom 
and the (— l)'dth cavity field. 

The total Hamiltonian H — Hj + He exihibits differ- 
ent behaviors in the long- wavelength (low-energy) regime 
and the short-wavelength (higher-energy) regime, which 
correspond to the quadratic and linear regimes of the 
photon dispersion relation, respectively. Namely, in the 
low-energy regime, the long-wavelength approximation 
gives a photon quadratic spectrum 



Ek~uj-2^ + ^k^ 



(6) 
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while in the higher-energy regime, the short-wavelength 
approximation leads to a photon linear spectrum 

Ek~uj-Tii± 2^k. (7) 

Both regimes will be studied in this paper. 

We note that an ideal system without losses is con- 
sidered here. In practice, both photons and atoms un- 
avoidably interact with different environments, that is, 
dissipation always exists. The dissipation substantially 
reduces the propagating length of the photons, and so 
does the transmission of the single-photon. In order to 
present the main physics of this system, wc neglect dis- 
sipation, decoherence, and the nonuniform couplings in 
this paper. These effects are separately studied and will 
be presented in the future. 

III. SINGLE-PHOTON REFLECTION AND 
TRANSMISSION 

A photon incident from the left of the CRW, with en- 
ergy within the energy band, propagates along the one- 
dimensional CRW until it is scattered by the first atom. 
Then it splits into a transmitted and a reflected portions. 
The transmitted part propagates freely until it encoun- 
ters the second atom, where the same type of splitting 
occurs once again. 

In this section, we will discuss the reflection and trans- 
mission coefficients of a single-photon in terms of the pro- 
jection of the asymptotic wavepackets onto appropriate 
plane waves. First, we consider the eigenstates of the 
total system. Three mutually-exclusive possibilities are 
considered: either the photon is propagating inside the 
cavity, or the photon is absorbed by one atom or the 
other. Considering all of three cases, the stationary state 
for the Hamiltonian H = Hi + He is written in the form 

\Ek) = E (j') 4 |0-95> + <e |Oe.g> + |Oge) , (8) 
j 

where the first number inside the Dirac brakets rep- 
resents the vacuum state of all cavity fields. The pa- 
rameter Uk {j) represents the probability amplitude for 
finding the photon at the jth cavity, ufj. is the probabil- 
ity amplitude of the Zth atom in its excited state while 
the other atom is in the ground state and all the cavity 
fields are in the vacuum. This form of \Ek) includes the 
three cases listed above. Using the Schrodinger equation, 
the single-photon scattering process can be described by 
the following equation 

{Ek - uj) Uk {]) + £, [uk U + l)+Uk {] - 1)] (9) 
= JiGik5j{-d)Uk i~d) + J2G2k5jdUk (d) , 

where the Green function Gik — Ji/ {Ek — fi;) is obtained 
from the relation 

< = ^;^^^>^[i-^fd\. (10) 



If we regard the second term of the left side of Eq. ([9]) as 
the kinetic energy term and the right side of Eq. ^ as 
potential energy term, then Eq. ([9]) describes the eigen- 
function Uk (j) subjected to a potential with singularities 
at j — ±d. In the region j ^ ±d, the potential is zero, the 
solutions to Eq. ^ are plane waves with wave-vectors k. 
Therefore we consider the wave functions 

r e^'^J + re-''^^ j < ~d, 
Uk (j) = <^ Ae'''^ + Be-'^^, -~ d < j < d, (11) 
[ te^^^ , j> d. 

The Uk (j) in Eq. ((TT|) describes the scattering process 
of an initial plane wave exp {ikj) incident from the left 
side of j = —d. This freely propagating wave is either 
refiected or transmitted when it encounters the first scat- 
terer. The reflection and transmission are described by 
the amplitudes r and A. The transmitted wave prop- 
agates freely until it encounters the second scatterer at 
the point j = d, where the corresponding reflection and 
transmission amplitudes are described by B and t. Here 
both scatterers produce a highly localized repulsive or 
attractive effective force, which depends on the incident 
energy of the single photon, the transition energies of the 
atoms, and the coupling strength between the atoms and 
their corresponding cavities. 



A. atomic transition energy ili inside the band 

Figured] schematically illustrates the seven cases (a-g) 
shown in Fig. [3l Figure [3] shows how the potential en- 
ergy depends on the energy Ek of the incident photon. 
All the results schematically shown in Fig. [3] are derived 
directly from Eq. Here, we assume fii < and 

both transition energies = 1,2) are in the region 
[uj — 2^,uj + 2^]. When the photon energy Ek < ^1 is 
smaller than the transition energy of the first atom, two 
attractive delta function potentials appear at j — ±d 
(as schematically shown in Figs, [^h- and [3]a,) . If the pho- 
ton energy Ek is increased and approaches fli, Figs.l^b) 
and[3]^b), the potential located at j = —d tends to mi- 
nus infinity. However, when Ek approaches fli from the 
right side, as shown in Figs, dfc) and El^c), the first po- 
tential gradually becomes an infinite barrier. When the 
incident photon energy Ek is further increased and is be- 
tween fti and fl2, Fig. [Ifd), photons first collide with a 
repulsive finite potential, then go through an attractive 
potential well, as schematically shown in Fig. ^d). As 
Ek further increases. Fig. [^Je), the height of the poten- 
tial barrier at j = —d becomes lower and lower, and the 
second two-level system creates a potential well with its 
depth becoming deeper and deeper, Fig.[3](e), eventually 
becoming a potential well with infinite depth. After Ek 
goes across the transition energy V,2, Fig.^f), a double- 
barrier is produced by the atoms, as shown in Fig. [2If)- 
In this case, the waves are totally reflected when the 
height of the delta potential of the second atom, goes 



(a) 



(c) 



(d) 



(e) i I (f) 



(g) 
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FIG. 2: (Color online). This figure schematically illustr 
seven different cases (a-g) discussed in the text and also 
responding to the seven cases (a-g) in Fig. [S] Here E 
the photon energy, and Q.i,Q.2 are the two atomic transi 
energies. The arrows indicate when the photon energy 
increases or decreases. In: (a) Ek < f^i; (b) Ek is slig 
below fii , and increasing; (c) Ek is slightly above fii , anc 
creasing; (d) €li < Ek < 0.2\ (e) Ek -^5^2; (f) Ek approa 
Q.2 from the right side, (g) Ek moving away from 572. 



to infinite, Fig.[3ff). Figure El^g) shows the potential 
ergy corresponding to the energies shown in Fig. [2] 
Therefore, Figure [3] schematically presents ways to ( 
trol photon transport by, e.g., adjusting the transi 
frequencies, fii and 0,2, of the two atoms. 

From the continuity conditions Uk {±d^) = (± 
and the eigenvalue equations 



{Ek-uj - J2Gk2) Uk [d) 
= {.d+ 1) + Uk{d^ 



1)], 



{Ek - uj - JiGki)uki-d) 

= -eK i-d + I) + Uk {-d - I)] 

at j ~ ±d, the transmission amplitude t can be derii 

t = 4e'sin2fc[(e^4^''^-l) JiGfeiJ2Gfc2 
+2i^ sin kJ^JiGki +4:f sin^ 

Above, d^ = d± e, where e is a very small positive num- 
ber. When the photon frequency Ek matches one of the 
atomic transition frequencies f2; = 1,2), the transmis- 
sion t is zero. This t — case occurs in cases (b), (c), (e), 
and (f ) in Figs. and [31 When an atom has its transi- 
tion energy fli inside the energy band, it may be excited 
by the incident photon. The absorption or emission of 
a photon by an atom leads to wave interference between 
the incident wave and the reflected wave. 



B. atomic transition energy Qi outside the band: 

0.1 < CO -2^,0.1 >Lu + 2^ 

For an atom with transition energy Qi far away from 
the energy band, the propagating single photon cannot 
excite the atom, thus the photon emerges in the other 
side of the scatterers with its energy (almost) equal to 
its original one, due to energy conservation. 

Using Eq. (|14p. in Fig. 21 we plot the reflection coef- 
ficient R = 1 — T (blue solid line) and the transmission 
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FIG. 3: (Color online). Schematic diagrams of the single- 
photon scattering process with Qi < ^2 inside the energy 
band. Here, the vertical axis is energy and the horizontal 
axis is the position along the CRW. Also, Qi is the transi- 
tion energy of the l-th atom (i = 1,2). The dashed lines 
refer to an infinitely high potential barrier or well. The seven 
cases (a-g) shown here correspond to the seven photon energy 
regimes (a-g) explained in Fig. [2] The cases (b), (c), (e), and 
(f) correspond to the total reflection of the incident photon 
because in each one of these cases one of the wells or barriers 
is infinite. 



coefficient T = |tp (red dashed line) as a function of 
the incident photon wave number k (— tt < k < tt), when 
both atomic transition energies f2; are outside the energy 
band. As shown in Fig. [4l total reflection R — 1 always 
happens at fc = 0, ±7r for nonzero Ji and J2, and this to- 
tal reflection is completely independent of the transition 
energies fli of the atoms. This observation is caused by 
the following reasons. First, a band, Eq. ([4]), is formed 
in this periodic CRW, which acts like a photon filter or 
a photon "band-pass filter" : transmitting photons over 
a limited frequency range. Incident photons with en- 
ergy Ek outside the band do not interact with the CRW, 
therefore the photon group velocities vanish at the zone 
boundaries or band edges. Second, when the atomic tran- 
sition energies f2/ are outside the band, the infinite delta 
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FIG. 4: (Color online). The photon reflection coeflicient R — 
1 — T (blue solid line) and the photon transmission coeflicient, 
T (red dashed line) as a function of the photon wave number 
k (— TT < fc < tt) when both atomic transition frequencies fli 
are outside the band [w — 2^,u! + 2^]. Here, k and d = 10 
is in units of the lattice constant, and other parameters are 
in units of ^. Here, the cavity energy is w = 5, (a) fii = 8, 
fl2 = 8, and the coupling between the atoms and the CRW 
are Ji = 0.5, J2 = 0.7; (b) Q.i ^2,0.2= 8, Ji = 0.7, J2 = 2; 
(c) fii =2,0.2= 8, Ji = 0.7, J2 = 2.6; (d) 0.^ =2,0.2= 2.7, 
Ji = 0.5, J2 = 3. Panel (a) corresponds to two potential wells 
at j = ±d. Panels (b) and (c) correspond to one potential 
barrier at j = — d and one potential well at j = d. Panel (d) 
corresponds to two potential barriers. 



potential wells or barriers cannot be formed, and the total 
photon reflection does not occur, except when fc = 0, ±n. 

The oscillations shown in T(fc) and R{k) in Fig.|4]orig- 
inate from the multiple interference of waves in the re- 
gion sandwiched by the two atoms. Comparing Figs.HI^a) 
andUfb), we find that as the couphng strengths J; be- 
tween the atoms and the CRW increase ( Ji from 0.5 to 
0.7 and J2 from 0.7 to 2), the oscillation in R{k) be- 
comes much larger when |fc| < 1, i.e., increasing the cou- 
pling strengths Ji magnifies the oscillations in R{k) and 



T(fc). Indeed the wave interference giving rise to the os- 
cillations in R{k) and T(fc) varies with five parameters: 
the energy Ek of the incident photon, the atomic ener- 
gies fli, and the couplings J; between the atoms and the 
CRW. The case shown in Fig.Ua) corresponds to two po- 
tential energy wells. The cases considered in Figs. IHb) 
and correspond to one barrier in j = —d and one 
well in j = d. The case considered in Fig. [HJd), there 
is a double-barrier. As is well known, the reflection and 
transmission coefficients are the same for a delta poten- 
tial barrier and a delta potential well, but their reflection 
and transmission amplitudes are different by a phase fac- 
tor. It is this phase difference that produces the clearly 
visible oscillations, due to interference shown in Fig. [H 

Figure 0] shows the complex dependence of i?(fc) and 
r(fc) as a function of the coupling strengths Ji, the hop- 
ping energy (or the inter-cavity coupling strength) ^, and 
the detunings, 61 — uj — fli, between the atoms and their 
corresponding cavities. If both coupling strengths J; 
are much smaller than the hopping energy ^, the hop- 
ping plays a leading role. In this case, and as shown in 
Fig. IUJa), the transmission T(fc) is quite large. When 
either Ji is larger than its corresponding detuning Si, 
and also larger than ^, the reflection dominates, as in 
Fig. mjd) . This result can also be found from the phase 
diagram (Fig. 4) in Ref. fl7|. Figures HJb) and[DJc) show 
the intermediate stage between nearly total transmission 
in Fig. m^a) and nearly total reflection in Fig. [IJd). In 
order to make this point somewhat explicit, we now ap- 
proximately write the transmission amplitude as 



sin fc 



JAkd i\ JlGl J2G2 



2^ 2^ 



(15) 



■ ; , -27 

sm fc > — — h sm fc 



sin^ fc 



2C 



1 



2^5i 2^62 



J2 

+i sin fc 7 — ^ + sin^ fc 



when the coupling strength Ji is larger than the hop- 
ping energy and smaller than the corresponding detun- 
ing (5/. Here, Gi = Ji/61. Equation p5|) shows that, 
when 3> 2^(5/, the reflection spectrum is much larger 
than the transmission spectrum, which coincides with the 
change shown numerically from Fig. [Hb) to Fig. IDJc). 
Notice that in the large detuning condition {Si ^ Ji), 
the magnitude Jf /Si is the shift of the energy levels due 
to the atom-cavity interaction (see Figs. [5] and [6]). There- 
fore, the relation between Jf /Si and the half-width 2^ of 
the band determines whether the reflection or the trans- 
mission plays a dominate role. This phenomenon will be- 
come much clearer from the energy-level diagrams shown 
in the next section. 

When the atoms embedded in ±d cavities are identi- 
cal, e.g., O = fii = U.2, J = Ji = J21 the transmission 
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coefficient T = |t| is derived from Eq. (|T4|) as 



wliere 



JGk sin {2kd) cos {2kd) 



2^2 sin^ k 



^ sin k 



(16) 



(17) 



We note that tlie maximum magnitude of tfie trans- 
mission coefficient T in Eq. can be acfiieved wfien 
one of the below conditions is satisfied: (1) the cou- 
phng strength J is much smaller than the detuning 
6ph = Ek — ^ between the energy E^ of the incident pho- 
ton and the transition energy Qi of each atom for a given 
nonzero ^; (2) the coupling strength J is much smaller 
than the hopping energy between adjacent cavities for a 
definite nonzero detuning Sph] (3) the condition 

2^ sin k 



tan (2fcd) 



JGk 



is satisfied. This third leads to a resonant tunneling ef- 
fect, which will be studied below. 



IV. SUPER-CAVITY ON RESONANT STATES 

The above results show that two atoms may act 
as a potential double-barrier. Any potential double- 
barrier ca,n produce a wavefunction localized in space (see 
e-g-, [13, mi)- Therefore, photons located in the range 
[— d, d], between the two barriers, may bounce back and 
forth. Thus, this double-barrier forms a resonator [13] ■ 
Photons can leak out of the resonator owing to the fi- 
nite width and height of the potential barriers. There- 
fore the localized state formed by this potential energy 
double-barrier is called a quasi-bound state or a resonant 
state [29, 30, 31]. A particle tunnels, through two energy 
barriers, when its energy matches (resonates with) the 
localized energy level. 

In the system we consider here, the photon propagating 
in this system encounters a double potential well/barrier, 
separated by a distance 2d. Such a potential can exhibit 
resonances. In this section, we derive the conditions for 
the photons to be trapped inside the resonator formed 
by the two atoms. Therefore, a quantum super-cavity 
can be formed by two atoms embedded in two separated 
cavities of the CRW. 



A. photon wave function in the super-cavity 

For simplicity, we now assume that the atoms at j = 
±d are identical, e.g. Ji = J2 = J and fli = ^2 — 
il. The eigenfunctions in this system satisfy the discrete 
scattering equation derived by Eqs. and ^ 



{Ek - uj) Uk (j) 



<[uk{j + l)+ukij-l)] (18) 
-JGk [S-j-dUk (j) + SjdUk (j)] ■ 



A resonant state is an eigenfunction of Eq. (fT8|) under 
the boundary condition that only outgoing waves appear 
outside the potential. Therefore we assume that Eq. pB]) 
has the following solutions 

Ce"*'^-' , j < —d, 

Uk (j) = { Abe^'^^ + ke^^''^, ^d<j<d, (19) 
De'''^, j>d, 

where the coefficients A/,, Bi,, C and D are the ampli- 
tudes for finding the particle in the state exp (zLikj) re- 
spectively. We also define the normalized amplitudes via 
the ratios 

. G ^ D 

'"^A' 

By imposing the continuity equation Uk (d^) ~ Uk (d ) 
and using the Schrodinger Eq. at the point j = d, 
we find 



h 

d- 



JGke''^'"^ 
2i^ sin k — JGk ' 

iS. sin k 
2i^sinfc — JGk 



(20a) 
(20b) 



Using the continuity equation Uk {—d'^) — Uk {—d ) and 
the equation (|18p at the point j — ~d, we have 



h = ('?!||^-l]e--- 

-',2kdi^sink 



JG 



k 



(21a) 
(21b) 



Obviously, Eq. (pPk ) and Eq. (|21ap must be equal, and 
the odd and even parities of the quasi-bound states in 
Eq. (fT9|) are included in 



JGk 



(22) 



where the plus sign gives the even parity and the minus 
sign has an odd parity. 



B. existence of quasi-bound states 

Let us define a parameter A = 2^/ J^. Now we solve the 
transcendent equation (|22p by a perturbation approach 
up to second order in parameter A, e.g. O (A^). First, 
let us assume that equation (122]) possesses a real solution 
for k only up to first order in A (This assumption will be 
proved later in this section). We set k = k^e for a real 
wave number. When fc is a real number 



cos {2ki-cd) = =F 1, 
wo; ^^ _u 2Csinfci.e 

Sm (ZKi-eCt) = ± 



JGk.. 



(23a) 
(23b) 
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Then the relation 

tan (2fcrcd) 



2^ sin kic 



(24) 



provides the condition for the existence of quasi-bound 
levels, which lead to the transmission coefficient T = 1. 
For a resonant state with an odd parity, the momentum 
k satisfies 



(25) 



due to the zero probability for finding the particle outside 
the barriers, where e is a small positive quantity and n 
is an integer. Substituting Eq. ([25]) into the right side of 
Eq. (f24|) and with the condition ^ 2^, the momentum 
of the resonant state can be approximately obtained as 



2d 



{5 - 2^ cos qn) sin q„ + O (A^ 



(26) 



where 5 — lo — Vt, and qn — ti-k j d. For even-parity states, 
the momenta of the resonant states are similar to that 
in Eq. ([26|) with n replaced by (n -1-1/2). Equation ([26]) 
implies that discrete levels appear in the energy band. 
In a similar way, we can derive the discrete energy of the 
resonant states with even parity. 

Although Eq. (|26|) gives the energy of a quasi-bound 
state, it fails to describe the behavior of the wavefunc- 
tion outside the sandwiched region. Indeed, a complex 
wave number fc must be considered in order to obtain the 
lifetime of a quasi-bound state. To do this, let us come 
back to Eq. (|22p . Here we show that, to second order 
in A, the imaginary part of k (which can represent the 
lifetime of the resonant state via its dispersion relation) , 
appears when the wave number in Eq. ()26|) is treated as 
a complex number. We obtain approximate analytical 
expression of the wave number 



k = qn~ ^Qn + dQl 



(27) 



[6 cos qn - 2^ COS (2g„)] + O (A^) , 



where 



Qn = -AS cos qn) sin g„. 
d 



(28) 



Thus, the lifetime of the quasi-bound state with wave 
number k is given by the imaginary part of Eq. (j27p . 



C. quantum super-cavity 

In Fig. [5l we show the numerically-obtained spatial 
distribution of the photon wave function along the CRW 
for a given k in Eq. ((27)) . In Fig. [6l we also give the cor- 
responding schematic explanation for the establishment 
of the super-cavity of Fig. O In Fig. [SI the eigenenergy 
Lo of the cavity at the site j = ±c? is shown by green 




FIG. 5: (Color online). The probability \uk{j)\'^ for find- 
ing the photon for a given k versus the position along the 
CRW. The parameters are in units of J and are set as fol- 
lows: d = 10,w = 10,n = 3,^ = 0.2, (a) = 10, (b) = 6, 
(c) Q, — 7. j is in units of the lattice constant. For a given n, 
q„ — rni/d. These were the inputs to equations (|27|) and ((28}, 
which provide and k. With this fc, Eqs. (|19H22p are used 
to obtain the Uk{j)'s shown in the Figures. 



dashed lines, and the atomic transition energy Q is shown 
by the red dashed lines. The eigenvalues of two dressed 
states are denoted by the symbols e±. The blue solid 
lines present the energy band formed by other resonators. 
From Fig. (Sj it can be found that a well-localized state 
appears in the sandwiched segment as long as the cou- 
pling strength J is much larger than the hopping energy 
^ and the detuning S = lo — fl, as shown in Fig. (HJa). 
In this situation, the coupling strength J plays a domi- 
nant role, and thus the coupling J shifts the energies at 
j = ±d to 



1 



(29) 



When uj — fl, the two strong J-couplings split the orig- 
inal degenerate energies of the cavity and the single 
atomic-excited state into two new dressed states at the 
point j = ±c?. These two dressed states are outside the 
energy band of the incident photon, as shown in Fig.(Sfa). 
Therefore, the photons will have a very low probability 
of going through the atoms, since the resonance condi- 
tion is not satisfied in Fig. (D^a). Thus, if a photon is 
initially located between the two atoms, it will remain 
there, bouncing back and forth from the atoms. The 
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CO S 8 

(a) 




I S £ 

(b) 

FIG. 6: (Color online). The corresponding schematic expla- 
nation for the establishment of the super-cavity in Fig. ([SJ. 
The eigenenergy u) of the cavities at j — ±d are shown by 
the green dashed lines, and the atomic transition energy Q 
is shown by the red dotted line. The eigenvalues of the two 
dressed states are denoted by the symbols e±. The blue solid 
lines show the energy band formed by the numerous other 
resonators. The arrows denote the shifting of the bare eigen- 
values towards their dressed values e±. Here, (a) tj = fi, (b) 
uj > Q. 

wavefunctions shown in Figs. [UJbjC) also indicate that a 
super-cavity can be formed, but the leakage of this quan- 
tum super-cavity is larger than in Fig.[5Ia). The reason 
for this large leakage in Fig. [5fc) and specially ^h) lies 
in the energy diagram of Fig. [H] Although the coupling 
strength is much larger than the hopping constant ^ in 
Figs. [5fb,c), the original eigen-frequencies, described by 
uj and f2 in Fig. EUb), are shifted in opposite directions 
to £-1-, but this shift amount (e+ — w) is still inside the 
band, therefore the tunneling process may appear with 
larger probability than the case in Fig. [Sfa). Comparing 
Fig. m^b) with Fig. [njc), it shows that the probability for 
a single-photon in the outside region in Fig.[5l[b) is larger 
than that in Fig. [5^c). Figures [5]Jb,c) further show the 
relation between the magnitude J^/S and the half- width 
2^. When 6 ^ J and J > the dominant photon-atom 
couplings approximately shifts the energy of the cavity 
to 

+ (30) 

Indeed, figures [n{b,c) show the change of the resonant 
states when J'^/S approaches 2^, e.g., the relation be- 
tween the dressed state e+ and the upper edge of the 
band. Obviously, the dressed energy level £+ is closer to 
the edge of the band in Fig.OJc) than the one in Fig.[ni[b), 
therefore, the probability is much smaller in Fig. EJc) for 




(a) (b) ^ 

FIG. 7: (Color online). Contour plots of the norm square of 
the wave function itfe(j) with odd parity, versus the position 
along the CRW and the atomic transition energies fi, where 
n is assumed to be inside the band. Here, £^ — 2, ui = 10, 
d = 8 for (a) and d = 12 for (b). 

a photon to be outside the sandwiched segment. There- 
fore, each atom plays the role of a partially-reflecting 
mirror. Here, we present a way to tune the leakage of 
the quantum super-cavity. 

D. quantum super-cavity with Q, inside the band. 

The super-cavity was studied above for large coupling 
strength J. As long as J is nonzero, a perfect reflec- 
tion appears when the transition energy is inside the 
band, namely, when the single-photon resonates with an 
atom [l7j . Therefore, a perfect super-cavity exists re- 
gardless of the magnitude of J. Of course, a perfect 
supercavity (r = 1) is an ideal limiting case. In reality, 
decoherence and losses will make the reflection coefficient 
r < 1. 

The photon trapping energy can be found analyti- 
cally, since Eq. holds exactly when the transition 
frequency f2 satisfies the condition 

n = LU -2^C08{qn/2) (31) 

and the corresponding resonant state has wave numbers 
k — qn/2. Thus, in this case, the two atoms form two 
mirrors with perfect reflection, which leads to a perfect 
super-cavity. Figure [7] shows the contour plots of the 
probability for a single-photon as a function of the co- 
ordinate j along the CRW and for the transition energy 
f2. Here, the atomic transition energy 17 is assumed to 
be inside the energy band of the CRW. Also only the 
wavefunction Uk{j) with odd-parity is depicted in Fig. [3 
which is the reason why the probability |ufe(j)p is zero 
at j = 0. If one regards the segment between the two 
atoms as a finite chain with N ~ 2d sites, the wave num- 
ber in this segment takes N discrete values, which give 
rise to the discrete energy levels. Figure [7] shows that 
when the transition energy f2 matches one of the dis- 
crete energy levels in the segment sandwiched by the two 
atoms, bound states appear and a quantum super-cavity 
is formed. These bound states are independent of the 
coupling strength J; however, a nonzero J is necessary. 
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E. quantum super-cavity made of superconducting 
qubits 

Using superconducting charge qubits (one kind of "ar- 
tificial atom") as an example, we now focus on the ques- 
tion on trapping and re-emitting photons in this unusual 
type of atomic resonator. It is well known that the tran- 
sition frequency of superconducting charge qubits can 
be controlled by both the voltage applied to the gate and 
the external flux through the SQUID loop [Ull Hi]. 
Let us assume that a photon with energy 

Ek = En^uj~2^coaqn (32) 

is initially in the (— d)th cavity. First, we tune the tran- 
sition frequency fl outside the energy band and consider 
a large detuning (fi — i?„). When the photon meets the 
first qubit, it passes the first qubit and moves freely be- 
yond the first qubit due to the large detuning. After 
the photon is inside the spatial range [— d, d] between 
the two qubits, the transition frequencies Qi = Q2 = ^ 
are adjusted inside the energy band, and satisfy Eq. ([3T|) . 
Therefore the photon would be totally trapped inside the 
super-cavity. We note that a tunable super-cavity could 
also be obtained by doping two A-type atoms inside the 
coupled-cavity array [35| . 

Based on the previous discussion in this paper, we 
can conclude the following: (1) the single-photon can 
be trapped in the region [—d, d] with a finite lifetime; 
(2) the single-photon can get out of the atomic resonator 
when the transition frequencies fl of these two atoms are 
not equal to the incident energy of the single-photon. 
In the appendix we will show that the existence of the 
photon bound states between the two atomic mirrors is 
independent of the magnitude of the transition energy 
il. We also conclude that a new cavity is formed by the 
two atoms separately embedded in the two cavities of the 
coupled-cavity array. Therefore, in analogy with super- 
lattices in solid state, we call this cavity a super-cavity 
and the atoms act as atomic mirrors. 



V. LONG- WAVELENGTH EFFECTIVE 
THEORY 

In this section, we show that the real part of the mo- 
menta of the quasi-bound levels in the low-energy re- 
gion can be obtained by expanding the sine and co- 
sine functions in Eq. (|26p as sinq„ « q„ and cosg„ ~ 
1 — 9^/2. Low-energy photons propagating along the 
resonator waveguide have long wavelengths. Under the 
long-wavelength approximation, a quadratic spectrum 

^Lu^+ ik^ (33) 

is found by expanding the cosine function around zero 
in Eq. ([7|), where the superscript L in E^ refers to the 
long- wavelength or lower energy regime and = w — 2^. 



By introducing the field operator 

/oc 
dk exp [ikx) Uk (34) 

with the commutation relation 

[^{x),^Hx')]^S{x-x'), (35) 
the Hamiltonian of the system in real space becomes 

H= dV(^«-ea^)^ + ^{17|e),(e| (36) 




where — \g)i (e| is the spin lowering operator of the 
Ith atom. Since the total number of excitations is con- 
served, we consider the storage of a single photon in the 
region separated a distance 2d by two 5-potentials. In 
the coordinate representation, the stationary state of the 
system 

/OO 
dxuk{x)^^{x)\Ogg) (37) 
-00 

W,T/|0eg)+u[™^|05e) 

is the superposition of a single photon (first term) and 
the single-excited states of the two atoms (second and 
third terms). The effective equation for the photon 

JGk [S {x - d) Uk (d) +S{x + d)uk i~d)] (38) 

= {^dl + E{:-co^)uk{x) 

is achieved from the eigenvalue equation H = 
\^k)- Two (5-potentials appear in Eq. ([55)1 along the 
direction of the photon propagation, one is located at 
X = —d and the other is located at x — d. The height 
of the potential is dependent on the energy carried by 
a single-photon. These two atoms divide the region of 
photon propagation into three zones: (I) x < —d; (II) 
—d<x<d; (III) X > d. The effective Hamiltonian 

Heff - - (dl (39) 

is valid in all three zones, and corresponds to free-particle 
Schrodinger equations, except for the replacement of Ej^ 
by {E^-co^) mH\E{;)=E^\Et). 

We now concentrate on the case E^ > uj^. From the 
standard boundary conditions that a wavefunction is al- 
ways continuous and its derivative is continuous except 
at points where the potential is infinite, we can derive the 
continuity equations for wave function Uk (x) in different 
zones, and the discontinuity of its derivatives (slopes) at 
the points x — ±d. According to the symmetry of the 
system, we assume that Eq. psp has the following solu- 
tion 

( S'le-*'''^, x<-d, 
Uk {x) = <^ e**^^ + Blc-''''', -d<x<d, (40) 
[ £'26^*=^, x>d. 
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FIG. 8: (Color online). The probability \uk{x) \ for finding the 
photon in space using the long-wavelength effective theory. 
The distance between qubits is 2d = 10. Other parameters 
are in units of J. ^ = 0.1, n — 3, u — 5 (a) ^1 = 3, (b) ^1 — 4, 
(c) Q — 7. For a given n, qn ~ These were the inputs to 
equations (|43|l and (|44p . which provide and k. With this 
k, Eqs. (|40H42|I are used to obtain the Uk{x)^s shown in the 
Figures. 



Using the same approach described in the section IV, we 
obtain the coefficients 



^1 = 
Bl = 



fc^e-*2fc<i ^ jqL gjj^ (2fcd) 



2k^ + iJG^ ' 
2ikS, - JGl ' 



(41a) 
(41b) 
(41c) 



with — J I [E]^ — il) and the condition for the exis- 
tence of the resonant states 



i2kd 



±^k {E^ -n)Ti. 



(42) 



Here the wave number k is complex. Under the condition 
Q = 2^2/ (J^d^) < 1, the wave number 



L\2 



k = qn-^Qn+diQ_ 

+i^Qn{Si + S(ql)+0{Q') (43) 



is approximately obtained, up to second order in the pa- 
rameter Q, for those states with odd parity, and 



Qi. 



(44) 



The super- index L in Eqs. ((43l - |44)) refers to the low- 
energy regime (long-wavelength approximation) studied 
in this section. The real part, Re(fc), and the imaginary 
part, Im(/c) in Eq. (|43p of a quasi-bound state provide 
the energy and the lifetime of this state via the disper- 
sion relation in Eq. ([55]) . It is clear that the Re(fc) can 
be obtained by expanding the sine and cosine functions 
around zero. Obviously, when k = qn and fl = lo^ + ^/c^, 
a perfect cavity is formed. In this case, the coefficients 
Si and S2 are zero. In Fig. [51 the probability for finding 
a photon in space is shown. As the transition energy 
varies, the effective potential induced by the two qubits 
changes from barriers to wells. It can be found that, as 
the depth or height of the delta potential becomes larger, 
the leakage of the super-cavity becomes smaller, which 
offers a way to control the leakage of the super-cavity 
by adjusting the energy level spacing of the two qubits. 
Therefore single-photons can be trapped. 



VI. SHORT- WAVELENGTH EFFECTIVE 
THEORY 

In the higher-energy regime, the short-wavelength ap- 
proximation leads to a linear spectrum Ek — lu^ + 2^\k\, 
with lUt^ = uj — TT^. Introducing the left (right) bosonic 
field operator f\{x) {(p^p^{x)), which creates a left- 
moving (right-moving) particle at x, the tight-binding 
Hamiltonian now becomes 



/oo 
dxLpl^ {x) Lpa {x) 



(45) 



a=R,L 



2i^ / dx 



(Pr {x) dxfR {x) - {x) dxipL {x) 

The left-moving and right-moving fields interact with 
these atoms respectively, therefore the interaction Hamil- 
tonian becomes 

/oo 
dxS {x + {-Yd) [^i {x) S[- + h.c] 

+f7^|e),(e|. (46) 

Although the Hamiltonian in the short-wavelength 
regime (linear dispersion regime) is significantly differ- 
ent from previous ones, the number of total excitations 
is still a conserved quantity. The stationary state for 
H = Hc + Hj with one particle excitation takes the form 

l^fc) = E / dxuka{x)pUx)\099) (47) 
+uie \Qeg) + U2e \Oge) , 
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where the first number in Dirac bracket represents the 
vacuum state of the cavity fields. Hereafter, the sub- 
index "S" in will refer to the short-wavelength approxi- 
mation regime, u^r {x) and UkL {x) represent the proba- 
bility amplitudes for finding the photon along the right- 
moving and left- moving direction at position x. More- 
over, Ujf, (with i — 1, 2) are the probability amplitudes 
for one qubit in the excited state and the other one in the 
ground state. From the Schrodinger equation, we obtain 
the relation between the left-moving amplitude and the 
atomic amplitude in the excited state 



(48) 

The relation between the right- moving amplitude and the 
atomic amplitude is 



{El- 



2i^dx) UkR ^ jy^^5^2 



(49) 

We can also find that the atomic amplitude Uje, the right- 
going amplitude Uku., and left-moving amplitudes UkL 
satisfy the relation 



dxS 



X + d {ukR + UkL) 



(50) 



with the Green function Gf = J/ (E'f — fl) . After ehm- 
inating the variables uie and U2e, both the left-moving 
eigenfunction 

{El ^Lo,- i2^d,) UkL (x) (51) 

/OO 
dx'6 {x' - d) [ukR (x') + UkL [x')] 
-OO 

/OO 
dx'6 {x' + d) [ukR {x') + UkL {x')] 
'OO 

and right-moving eigenfunction 

{El -L0^ + i2^d,) UkR (x) (52) 

/OO 
dx'6 {x' - d) [ukR {x') + UkL (x')] 
-OO 
pOO 

+JGi6{x + d) / dx'6{x' + d)[ukR{x') + UkL{x')] 



are subjected to a delta potential with singularities at 
X — zLd. 

In the region x ^ ±d, the potential is zero, and 
the solutions of Eqs. ((5T|) and ((52)) are plane waves 
with left-moving and right-moving wave-vector number 
k = E^ /vg. Therefore, we can assume the right-moving 



UkR (x) 




(53) 



and the left-moving wave-function 



UkL {x) 



tL e 
tl e" 
0, 



-ikx 



X < — d, 
d < X < d, 
X > d, 



(54) 



which allow the existence of quasi-bound states in this 
system. The magnitude of tl 



■Tl 



JGk ^i2kd 



i2^ - JGl 
and the relations 

tR 



- JGt^_,^kd 



tle 
rL + e 



i2kd 



-i2kd 



(55) 



(56a) 
(56b) 



of the amplitudes Ir^ II and tl can be obtained by 
integrating Eqs. ([STjl and (|52|) in the neighborhood of 
X = zLd. For the appearance of quasi-bound states in the 
spatial range sandwiched by two atoms, Eq. ([55]) leads to 
the condition 



^2ikd _j_ 



2i^ 



(57) 



with the complex wave number k. Here, the lower sign 
corresponds to the odd-parity, and the upper sign corre- 
sponds to the even-parity. Obviously, when the transition 
energies of the two atoms are 



2^ 



(58) 



the bound states have odd parity. However the even par- 
ity corresponds to the transition energy 



n = 



2i 









4) 







(59) 



Except the situation discussed above, Eq. ([57|) does not 
have an exact solution. We now seek the values of k 
for which Eq. ([57)1 can be approximately solved. Here 
we only consider the energy levels with odd parity. A 
similar calculation provides results for even parity. Using 
the approach described above with the parameter P = 
4^2/(rfj2) and \ = 2^/j2, Eq. ^ with the lower sign, 
yields the wave number 



k 



J 

2d^ 



'd{Qfy+i^^Q^^+0{P') (60) 

whose real part can be obtained from Eq. ((26|) by ex- 
panding the sine and cosine functions as sinq„ sa 1 and 
cosQn around tt/2. Here, 



A 



{6^ + 2^qn) 



and 



6. 



7r£,-n. 



(61) 



(62) 



We plot the norm square of the left-going wave-function 
in Fig. IDJa), the right-going wave- function in Fig. (H^b), 
and the total wave-function Uk{j) in Fig. [IJc), where 

Uk{j) = UkL{j) + UkRij). 
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(b) 



(0 



FIG. 9: (Color online). The norm square of the left-moving 
wavefunction \ukL{j)\'^ (a), the right-moving wavefunction 
WkR{j)f (b), and probability |ttfe(j)P = \ukL{j) 
(c) for finding the photon in space. The parameters are set as 
follows: d = 8, f = 0.1, n = 1, u — 5, ^ — 2. Parameters are 
in units of J. For a given n, q„ = ^ . These were the inputs 
to equations (I60[l and (|61[). which provide and k. With 
this k, Eqs. (|53H57|I are used to obtain the iife(a;)'s shown in 
the Figs.(a)-(c) 



VII. CONCLUSIONS 

We have studied the coherent control of single pho- 
ton transfer in a coupled resonator waveguide with two 
atoms. The coherent control can be realized by adjust- 
ing the detuning between the single photon frequency and 
the energy-level-spacings of the atoms. We have shown 
that a super-cavity is formed in the coupled-cavity ar- 
ray due to the strong coupling between the atoms and 
the corresponding cavities, and the discrete values of the 
photon momenta are analytically derived. Moreover, a 
perfect super-cavity appears when the transition ener- 
gies of the two atoms are equal to the energy of an inci- 
dent photon. We also find that besides the bound states 
formed by two perfect atomic mirrors, there always exist 
other bound states at the edge of the band. The real 
parts of the discrete momenta obtained by the discrete 
approach unify those obtained by the effective continuum 
theory in both the long- wavelength and short-wavelength 
regions. 
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APPENDIX A: PHOTON BOUND STATES 
BETWEEN TWO ATOMIC MIRRORS 



To find the wave-functions for the eigenvalue equa- 
tion (jlSp . one needs to write down the wave-functions 
in different regions. Since exchanging of the two atoms 
does not change the equations for the wave-functions of 
the photon propagating along the CRW, here we only 
deal with odd-parity wave functions, which have the sinh 
function in the center region and exponential decay in 
the edge regions 

-A exp [{inn + k) j] x < —d, 

(.t) = ^ B exp (imrj) sinh (kj) — d < x < d, 
A exp [{imr — k) j] x > d. 



(Al) 



From the continuity and discontinuity conditions at x 
d, 

u u (d-) , 

+ JG^ - E)u{d) ^ i[u{d + I) + u{d - 1)] , 



we can easily obtain 



tanh {nd) 



^ exp (—imr) sinh k 



E-Lu- JG + ^{e 



cosh k) 
(A2) 



with 



and 



J 



E^-n 

In principle, n can be obtained by solving the implicit 
transcendental equation (|A2[) . It is obvious that k = is 
one of solution of Eq. (|A2|1 . This k = solution makes 
sure that the odd-parity wave functions exist and two 
bound states appear at the edges. 
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